Relaxation mechanism for electron spin in the impurity band of n-doped 

semiconductors 
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We propose a mechanism to describe spin relaxation in n-doped III-V semiconductors close to 
the Mott metal-insulator transition. Taking into account the spin-orbit interaction induced spin 
admixture in the hydrogenic donor states, we build a tight-binding model for the spin-dependent 
impurity band. Since the hopping amplitudes with spin flip are considerably smaller than the spin 
conserving counterparts, the resulting spin lifetime is very large. We estimate the spin lifetime from 
the diffusive accumulation of spin rotations associated with the electron hopping. Our result is 
larger but of the same order of magnitude than the experimental value. Therefore the proposed 
mechanism has to be included when describing spin relaxation in the impurity band. 

PACS numbers: 72.25.Rb, 71.70.Ej, 71.30.+h, 71.55.Eq, 



The renewed interest in spin relaxation in 
semiconductors 1 - stems from the possible applica- 
tions and fundamental science associated with the 
emerging field of spintronics 2,3 . The large measured 
values of the electron spin lifetime constitute a promise 
for the use of spin as a unit of quantum information and 
pose a considerable challenge for the identification of 
the appropriate mechanisms of spin relaxation. Interest- 
ingly, low-temperature experiments in various n-doped 
semiconductor bulk systems have found that the longest 
spin lifetimes occur close to the Mott metal-insulator 
transition (MIT) densitj*«i2jl. Many aspects of the 
Mott transition have been thoroughly studied, making 
it a paradigm of Condensed Matter Physics^. However, 
the connection between the spin and transport problems 
is only beginning to be explored 9 . 

The comprehensive experimental and theoretical work 
of Dzhioev et al. 7 considered both sides of the MIT, which 
for GaAs occurs at n c ~ 2 x 10 16 cm -3 . In the deeply 
localized regime, with donor densities n < 5 x 10 15 cm" 3 , 
the hyperfine interaction was shown to account for the 
measured spin lifetimes^. For higher densities, but 
still lower than n c , the anisotropic exchange of local- 
ized spins was proposed as the dominant mechanism for 
spin relaxation^. Later calculations based on the same 
mechanism- 2 initiated an ongoing controversy concern- 
ing the quantitative agreement between experiment and 
theory. 

For densities above the transition, the well-known 
D'yakonov-Perel (DP) mechanism for conduction elec- 
trons was invoked^. This mechanism arises from the 
splitting of the conduction band due to spin-orbit inter- 
action and yields a spin lifetime inversely proportional 
to the momentum relaxation time. This description ap- 
plies to doping densities high enough so that mainly the 
conduction band is populated. (For GaAs the hybridiza- 
tion of impurity and conduction bands occurs at a doping 
density n h ~8x 10 16 cm" 3 fiTH^.) If one is interested 
in understanding the large lifetimes measured near the 



MIT, the DP mechanism is not applicable, since, clearly, 
any mechanism invoking momentum relaxation via im- 
purity scattering becomes meaningless in the impurity 
band. This difficulty lies at the origin of the lack of suit- 
able theories of spin relaxation at low temperature near 
the MITi£. 

In this work we provide a theoretical frame and propose 
a spin-relaxation mechanism for the metallic-conduction 
regime of the impurity band. Our approach is to ex- 
tend, by incorporating the spin-orbit interaction, the 
well-known model of Matsubara and Toyozawa (MT) of 
electron conduction at zero temperature^. The effect of 
the spin-orbit interaction is to introduce spin admixture 
in the impurity states. We will refer to this as impurity 
spin admixture (ISA). The ISA allows for spin-flip pro- 
cesses in electron-hopping events even in the absence of 
spin-dependent potentials. A tight-binding model built 
on the ISA states provides a theoretical framework to 
study spin dynamics and spin-dependent transport in the 
impurity band. Since we are interested in spin lifetimes, 
a first test of our model is to estimate the order of mag- 
nitude that the ISA mechanism predicts. We proceed by 
calculating the accumulated spin rotation angle along the 
diffusive evolution of the electron in the potential of the 
impurities. The time that it takes the spin to depart an 
angle of unity from its initial orientation is then taken as 
a qualitative measure of the spin lifetime. 

The MT model consists of a tight-binding approxima- 
tion built from the ground state of the doping impuri- 
ties. For shallow donors it is a standard approximation 
to restrict the expansion of the impurity ground state to 
conduction-band stated and for an impurity located at 
the origin we write it as 



hM(r) = 5>(k)e*' r [u kl<r ](r) « 0(r) 
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The envelope function tfi(r) = (l/ira 3 ) 1 ' 2 exp (— r/a), 
where a is the effective Bohr radius (a « lOOA for GaAs), 
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is hydrogen-like. We note as </>(k) its Fourier transform, 
while [uk,a] ( r ) represents the periodic part of the Bloch 
functions of the conduction band states. Its dependence 
on k, being much smoother than that of <fi(k) , leads to the 
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last relation in Eq. (p}, where we note 

The spinors [ipg, CT ] and [lik.o-] are trivial since they are 
eigenstates of the operator S z with eigenvalue a = ±1. 
However, this will no longer be the case once we include 
the spin-orbit interaction. The Hamiltonian of the MT 
model can be simply written as 
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where c^,^ (c m(T ) represents the creation (annihilation) 
of an impurity eigenstate at the impurity site ml (m). 
The ground-state energy of an isolated impurity is taken 
as the origin of energies. The energy integral for the 
electronic transfer from site m to ml is given by 
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The Coulomb-like potential produced by the impurity 
placed at r m is V m (r) = —e 2 /e\r — r m |. We note e the 
static dielectric constant (12.9 for GaAs), Vo — e 2 /ea, e 
the electron charge, and r mm i the distance between the 
two impurities. For convenience, in Eq. ([3]) we switched 
from spinor to ket notation. The Hamiltonian Eq. @ has 
been thoroughly studied using a variety of analytical and 
numerical technique s 14 ' 18 ! 19 ' 20 , allowing a useful descrip- 
tion of the impurity band and its electronic transport. 

In order to extend the MT model to the spin case 
we first generalize the shallow-donor wave functions to 
include the spin-orbit interaction. The spin-orbit ef- 
fects coming from the orbital motion do not modify in 
an appreciable way the envelope functions 4>(r). There- 
fore, the spin-orbit interaction affects mainly the spinor 
[uk]- As is well-known, in lattices lacking inversion sym- 
metry (i.e. zincblende semiconductors like GaAs), the 
spin-orbit coupling leads to spin-mixed conduction-band 
states. Group-theoretical arguments dictate the way in 
which the conduction and valence states are mixed by the 
spin-orbit interaction. Within the k ■ p approximation of 
Kane 21 (and using the notation of Ref. 22), the periodic 
part of the spin-mixed conduction-band states is given 
by 



«k CT ) = k 0) )+k-|uy 



where 



uW) = a x (|Ro) + a 2 S x |Ro)) 
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The state \u a } is s-like, since it describes the un- 
perturbed wave function at the T-point. The vector 
R) = (\X), \Y), \Z)) represents the three p-like valence 
states. S is the spin operator. The state |uko-} is clearly 



not an eigenstate of S z . However, we still character- 
ize it with the label a since the mixing is small, and 
{uk<j\S z \uka) is much closer to ah/2 than to —ah/2. 
The weak spin mixing is governed by the small con- 
stants ax = ih[(3E G + 2A)/{6m*E G {E G + A))] 1/2 and 
a 2 = 2A/ih(2A + 3E G ). We note A the spin-orbit split- 
ting of the valence bands, m* the conduction-band effec- 
tive mass, and E G the bandgap. 

The mixing of Eq. tells us that in the presence of 
SO interaction, Eq. (JTJ) takes the form 
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Using the hydrogenic character of 4>(r), the impurity 
spin- admixture (ISA) state centered at r m reads 
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Electron hopping between ISA states in different impu- 
rity sites provides a mechanism for spin-flip by connecting 
the a and er = — a states. The Hamiltonian of the system 
now contains a term without spin-flip (described by the 
Hamiltonian Hq of Eq. (f5|)) and a spin- flip term 
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As in the spinless case, the matrix element 
(ipm'crlVrn'lipmo-) is expected to dominate the energy 
integral for the electronic transfer between states m and 
ml . With spin flip, however, symmetry reasons dictate 
that the corresponding matrix element {ip m 'a\Vm' IVVncr) 
vanishes. This important fact is ultimately responsible 
for the large values of the spin lifetime in the regime 
of impurity-band conduction. We are then forced to 
consider the matrix element 



(.WIWrrKr) = ~ C I ^ ( 9 ) 

(r - r m )-(z - z m ,) - (z - z m )(r - r m /)_ 



r - r m r - r„ r - r m / 
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with p m, ml, r± = x ± iy, C = Vo\ai\ 2 a3/ira 4 , and 
a 3 = 3A(A + 2£ G )/(2A + 3E G ) 2 . Three-center inte- 
grals like the one of Eq. ([9]) are in general very difficult 
to calculate 2 ^. However, in the case r mm >/a 3> 1 H3 we 
can perform the integrals using the saddle-point approx- 
imation and obtain 

(^m'al^lVw) = -4.2 Ce lVm (cos0 p + i cos9 m sin^ p ) 
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ip m and 6 m are the polar angles of the vector r mm i in 
the original coordinate system. <j) p , p p , and z p are the 
cylindrical coordinates of r p in a new system having the 
z-axis parallel to r mm / and its origin at the middle point 
between the sites m and m! . For a given impurity con- 
figuration, the matrix elements (|10|) yield the energy in- 
tegral with spin flip 
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of the Hamiltonian Hi (Eq. ©). The system Hamilto- 
nian Hq + Hi can be addressed numerically or by pertur- 
bation theory. In order to test the physical relevance of 
the proposed spin-flip mechanism, we will estimate the 
spin relaxation time within some simplifying hypothesis 
that we discuss in what follows. 

Viewing the electron transport as a hop between im- 
purity sites, we see that, since \V^ m , \ ^> \tmm'\> there is 
a very small probability of spin-flip per hop, which may 
be translated into a mean spin-rotation angle 7 m m' . As- 
suming that these relative rotations are accumulated in 
a diffusive way, we can estimate the spin-relaxation time 
by 
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where < 7 2 > is the ensemble average of 7 mm / and r c is 
the mean hopping time. If in the hop between impurities 
to and m! the electron is initially in the spin-up state, 
the expectation value of the spin operator after the hop 
is 
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The angle between the initial and final spin orientations 
is 



v = Arccos 
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Allowing for an arbitrary initial spin orientation be- 
fore the hop enhances the root mean square value of 



Omm' by a factor of 
(3/2) 2 < 6 



'2. We then have < 7 mm , >= 
>, where an additional factor of 3/2 ap- 
pears since 9 m m' contains the two components of 7mm' 
that are relevant for spin relaxation. 

Given the form of (TIT)) of the energy integrals, the typ- 
ical rotation angle 7 mm / involves a double sum over im- 
purities p and p' (7^ to, to'), which can be approximated 
by its impurity average. Only the diagonal term (p = p') 
survives the average, yielding 



7mm' — 33.2 ( ; ; 



3/2 11/2 
'mm' u 



ni(l + cos 2 6> mm /), (15) 



where n\ is the impurity density. The typical rotation 
angle between impurities to and to' increases with their 



distance r mm / as a power law. This dependence renders 
the impurity distribution crucial for the determination of 
the mean square rotation angle per hop. The distribution 
of doping impurities is known to be completely random 
and to lack hard-core repulsive effects on the scale of 
a [251 ]. Since the probability of jumping from a given 
impurity to to a second one m' is iCm'lVEm' l*wl 2 > 
we obtain the impurity average of the typical rotation 
angle per hop as 



< 7 2 >= 1.8 x 10 2 ( — 



a (n;a 
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The hopping time can be estimated from perturbation 
theory by determining the characteristic time for the de- 
cay of the initial population by one-half, yielding 
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(17) 



where an average over the initial position to is implicit. 
This estimation assumes orthogonality of the electronic 
orbitals at different impurity sites, which is actually not 
satisfied by the hydrogenic states. However, since the 
overlaps are very small, the non-orthogonality effects 
arising in the MT model are known to be smal l 18 ' 26 . 

We calculate the rotation angle for hopping processes 
between impurity sites. However, if we take delocalized 
impurity-band initial states, a golden rule approach with 
some simplifying assumptions leads to the same expres- 
sion for r s up to a factor of order one2£. 

From Eqs. |T2]), HI]), and {TTJ) we obtain 
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For GaAs, at the density of the MIT, Eq. (18]) yields a 
spin-relaxation time of r s = 1200 ns. 

This value is larger but within an order of magnitude 
of the experimentally reported result of 200 ns 6.7] at the 
MIT. We stress that our result does not depend on any 
adjustable parameter, but it relies on a few approxima- 
tions. For example, a step of the calculation that could 
admit an alternative treatment is the impurity average of 
7mm' • 11 we assume that the hopping takes place only be- 
tween nearest neighbors (separated by a typical distance 

— 1/3 

n i ) we obtain a value of 700 ns for t s at the MIT. 
In general, the order of magnitude agreement with the 
experimental value is not affected by the approximations 
and therefore the ISA mechanism needs to be included in 
the description of spin relaxation in the impurity band. 

In materials with stronger spin-orbit interaction like 
InSb and InAs, one obtains considerably smaller val- 
ues of the spin lifetime. For the impurity density of 
5 x 10 14 cm~ 3 in InSb, our estimate yields a spin life- 
time of 86 ns. This value is within an order of magnitude 
of the experimental result and close to the theoretical 
estimation of Ref. [22|. The theoretical approach of this 
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reference is justified only in the high-concentration limit 
since it describes single scattering of plane waves and 
treats electron-electron interactions through a corrective 
factor. 

In this work we have proposed a new mechanism for 
spin relaxation in the regime where electron conduction 
occurs in the impurity band of doped semiconductors. 
The mechanism is based on the impurity spin admixture 
of the electronic ground state of the donors caused by 
the spin-orbit interaction. The impurity spin admixture 
states do not have a well-defined spin projection about 
a fixed spatial direction, and therefore hopping between 
two of these states may connect different projections of 
the angular momentum. 

Unlike the spinless case, the matrix elements of the 
spin-flip hops are not dominated by a two-center integral, 
where the impurity potential corresponds to one of the 
extreme sites. We therefore have to consider three-center 
integrals, where the impurity potential is not centered 
around any of the two sites of the hop. Since the latter 
matrix elements are considerably reduced with respect to 
the former, the resulting spin lifetimes are very large. 

Our calculation of spin-flip matrix elements yields a 
suitable model for studying electron and spin transfer in 



the regime of impurity concentrations just above that of 
the metal-insulator transition. Various treatments can 
be applied to our model Hamiltonian. In this work we 
provide an estimation of the spin lifetime by calculating 
the diffusive accumulation of spin rotation during the 
hopping process. This estimation yields values that are 
larger than the ones experimentally measured but within 
the right order of magnitude. Therefore the impurity spin 
admixture mechanism has to be taken into account in 
descriptions of the spin relaxation in the impurity band. 
Our model admits generalizations including other phys- 
ical effects, like doping compensation, electron-electron 
interaction, and a second electronic band, which may im- 
prove the agreement between theory and experiment. 
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